Filtered topological structure of the QCD vacuum: effects of dynamical quarks 



Falk Bruckmann, Florian Gruber, Andreas Schafer 

Institut fur Theoretische Physik, Universitdt Regensburg,D-93040 Regensburg, Germany 



Abstract 



o 

(N 



We systematically compare filtering methods used to extract topological structures on SU(3) lattice configurations. We 
show that there is a strong correlation of the topological charge densities obtained by APE and Stout smearing. To get 
rid of artifacts of these methods, we analyse structures that are also seen by Laplace filtering and indeed identify artifacts 
for strong smearing. The topological charge density in this combined analysis is more fragmented in the presence of 
dynamical quarks. A power law exponent that characterises the distribution of filtered topological clusters turns out to 
be not far off the values of an instanton gas model. 
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1. Introduction 

The distribution of topological charge density is a char- 
acteristic signal for topological objects and their appear- 
ance in the QCD vacuum. It is interesting to study also 
because the topological susceptibility is related to the rf- 
mass. 

Such non-perturbative aspects of QCD are best stud- 
ied through lattice simulations. However, the topological 
density in those configurations is - like many other observ- 
ables - dominated by quantum fluctuations at the scale of 
the lattice spacing. Many methods have been developed to 
extract the IR content from lattice configurations. Among 
them are cooling/smearing, Laplace filtering (for details 
see below) and the fermionic definition of the topological 
charge via eigenmodes of a chiral Dirac operator [TJ [2]. 
Unfortunately, all these methods introduce their own am- 
biguities and parameters, which may lead to wrong con- 
clusions. Cooling for instance is biased towards classical 
solutions after many steps. 

Hence, in order to get a coherent picture of the topolog- 
ical structure of the QCD vacuum, it is necessary to con- 
trol or even remove these ambiguities. A way to achieve 
this is the comparison of different filtering methods. In an 
earlier work [3J we have demonstrated that APE smearing, 
Laplace filtering and Dirac filtering can indeed be matched 
to each other on quenched SU(2) configurations. More- 
over, they agree very well on the 'hot spots' in the filtered 
topological charge, which in addition are subject to an in- 
teresting power-law. A similar comparison using overlap 
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(valence) quarks has been performed in Ref . [J] . Concern- 
ing the total topological charge from different definitions, 
the agreement has been shown to become almost perfect 
in simulations with actions that generate smoother config- 
urations a priori (within the Monte Carlo chain) thereby 
fixing the topological sector [5j [3] . 

In the present Letter we apply the methods of Ref. [3J 
to 517(3) configurations (which can be done with slight 
modifications) with the aim to study the difference of topo- 
logical structures in quenched vs. dynamical configura- 
tions. One might expect that light quarks suppress topo- 
logical objects due to their topological (near) zero modes. 
On the other hand, the Adelaide group has found an in- 
creased density of smaller nontrivial field configurations 
after fitting the filtered topological density to instanton 
profiles [TJ. This study, however, uses long smearing with- 
out comparison to other filtering methods and staggered 
quarks, whose topological/chiral properties are question- 
able. 

We will show that APE and Stout smearing yield very 
similar topological charge densities, i.e. agree on a non- 
perturbative observable, when the parameters are matched 
as in perturbative studies [HI H] • Nevertheless, to control 
the ambiguities it is necessary to complement smearing 
by an independent filtering method, in our case Laplace 
filtering. This comparison clearly reveals that smearing 
has to be used with care not to destroy small topological 
objects. 

Our main conclusion concerning the effect of dynamical 
quarks is that they lead to a more fragmented structure of 
the QCD vacuum. This is based on the power-law coeffi- 
cient and the absolute number of topological clusters com- 
mon to both methods. Therefore, this statement is free 
of ambiguities and also does not use a particular model of 
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topological objects in the QCD vacuum. 

The Letter is organized as follows. In the next two sec- 
tions we briefly review the definitions of the filtering meth- 
ods and the topological density we use. Sect. [3] deals with 
the comparison of the filtering methods and how they can 
be matched optimally. In Sect. [5] we analyze the filtered 
topological clusters and reveal differences in quenched vs. 
dynamical configurations. Finally, we give our conclusions. 

2. Filtering methods 

One of the first attempts to filter out the UV 'noise' 
has been APE smearing [101 [H] , defined as: 
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where a determines the weight of the old link and the 
sum of the attached staples. Throughout this work we 
have used 4-dimensional APE smearing with the standard 
value a = 0.45. 

The right-hand side has to be projected back to the 
gauge group. The projection onto SU (3) is not unique, but 
a common choice is to define Psu(3)(W) as the that ele- 
ment of the group V € SU(3) that maximizes Re Tr{VT0}. 
In our work, the maximum is found iteratively using code 
included in the CHROMA software package for lattice 
QCD P2]. 

Stout smearing [TJ] avoids this projection by using the 
exponential map: 
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is a hcrmitian traceless matrix constructed from all plaque- 
ttes containing the old link C/° ld and a smearing parameter 
As e l( 2 <E SU(3) we do not need any projection onto 
the gauge group. The identification of the hermitian part 
of a matrix ft is not unique, we use k(Cfl — fi). Both 
smearing methods are iterative procedures, which lead to 
smoother configurations in each step. These methods had 
been used with great success in previous studies of the 
QCD vacuum pj Q31 [IS]. But it is hard to decide how 



1 Here we also refer to a simultaneous smearing of spatial and 
temporal links like in the case of APE smearing 



many smearing steps one can apply without loosing infor- 
mation about the underlying structure. 

A relatively new and substantially different method is 
Laplace filtering |16j . This method is based on a spectral 
decomposition of the links in terms of cigenmodes $ n {x) 
of the covariant lattice Laplacian: 



Ntat 

n=l 



(4) 



where the sum runs over all Ntat = Yol-N c eigenmodes. As 
the eigenvalues are related to some energy scale (squared) , 
one can truncate the sum at N -C N tot to get rid of the 
high energy contributions: 



U^\x) = P su{3) 
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This procedure acts as a low-pass filter in the sense of a 
Fourier decomposition: with lower values TV of included 
modes the Laplace filtering gets stronger. Taking the fil- 
tered links as a starting point, one can in principal measure 
any observable. 

While both strategies aim at smoothing the gauge fields, 
it should be stressed that Laplace filtering is completely 
different from smearing, because it is based on extended 
objects, namely the eigenmodes, and does not locally mod- 
ify the gauge links in contrast to smearing. 

In this sense this method is quite similar to Dirac fil- 
tering [2], which provides direct access to the topological 
charge density via the low-lying eigenmodes of the Dirac 
operator. However, Laplacian modes are chirality blind 
and do not obey an index theorem, which connects the 
zero-modes of the (chiral) Dirac operator to the topolog- 
ical charge. From a computational point of view Laplace 
filtering is much cheaper than Dirac filtering (with good 
chirality properties). 

3. Topological Charge Density 

Throughout this Letter we will focus on the topological 
charge density: 



(6) 



where F pv is an improved discretization of the field strength 
tensor [17], which combines lxl, 2x2 and 3x3 loops to 
achieve 0( a 4 ) -improvement at tree-leveQ This definition 
is better behaved in the continuum limit than the naive 
discretization with only lxl loops and gives topological 
charges Q = J2 X q(x) closer to integers. 

In order to get results for dynamical and quenched con- 
figurations which can be easily compared, we have chosen 
two SU(3) ensembles with the same lattice spacing and 
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lat. 


size 


lat. spacing 


Plw 


quenched 


16 3 


• 32 


0.148 


7.90 


dynamical 


16 3 


• 32 


0.150 


4.65 



Table 1: Details on the gauge configurations used in this work. 




(a) APE smearing 




(b) Laplace filtering 



Figure 1: Comparison of the same slice of the topological charge 
density resulting from 7 APE smearing steps (a) and 500 Laplace 
modes (b). 

the same physical volume (see Tab. [I]). The ensembles 
were generated with the Liischer-Weisz gauge action and 
a chirally improved Dirac operator [TH]. For the dynami- 
cal simulations two flavors of mass degenerate light quarks 
were used (details can be found in [T9]). 

4. Comparison of the Methods 

For quenched SU{2) lattice configurations we had ob- 
served a similarity of the topological density obtained with 
the different filtering methods [3] . A corresponding exam- 
ple of a quenched SU(3) configuration after APE smearing 
and Laplace filtering is shown in Fig. [I] 

For a quantitative comparison one needs a measure of 
the similarity. In [3] the following quantity was introduced: 



XAA XBB 

with 

xab = (i/v) M x ) - M*) - 5fl) (8) 

x 

being the correlator of two topological charge densities 
<Ia{ x ) an d Qb(x) (mean values q = Q/V are subtracted 
for convenience). 

H is a positive quantity and equals 1, if the densities 
differ only by a constant scaling factor. The great advan- 
tage of this quantity is that we do not have to know the 



normalization factors of the topological charge densities, 
since these factors drop out in Eq. ^ and we are able to 
compare only the relative difference at each lattice point. 

For a quantitative analysis of the local structure result- 
ing from different methods A and B, it is useful to define 
'best matching' pairs of filter parameters, by maximizing 

S-AB- 

To be precise, we keep the smearing parameters con- 
stant and compute Sab for up to 50 APE steps and 50 
Stout steps respectively up to 500 Laplace modes. More 
smearing steps have not been applied, as too much smear- 
ing will destroy the local topological structure, whereas 
more Laplace modes would be too expensive. 

We first discuss the matching of APE and Stout smear- 
ing, see Fig.[2]upper row. In the absence of smearing (lower 
left corner) one has H = 1 as the configurations are iden- 
tical. 

An interesting observation is that the maximal values 
of S lie on some kind of ridge line, where 5 > 0.95 within 
our range of smearing steps. Thus, the different smearing 
techniques reveal almost the same structures. The slope 
of this line is determined by the ratio of the smearing pa- 
rameter. 

An almost one-to-one correspondence for the topolog- 
ical density after the same number of APE and Stout 
smearing steps (i.e. a slope of 1) is achieved for the pa- 
rameters related as a — 6 • p [20] . 

This is consistent with the perturbative result of Ref . [5] 
and with the non-perturbative result from Capitani et al. 
[5]. The latter group has focused on rather global observ- 
ables, which are not related to topology, with up to three 
smearing steps. Our non-perturbative result reflects the 
local similarity of both methods and their strongly corre- 
lated topological charge densities up to 50 steps. 

The described equivalence applies to both dynamical 
and quenched configurations. 

The comparison is not so perfect when we consider 
APE smearing and Laplace filtering, see Fig. [2] lower row. 
Overall, these plots are similar to Fig. 3 (bottom panel) 
of Ref. [3J. Although there is a qualitative similarity for 
quenched and dynamical configurations, we find a quanti- 
tative difference. The agreement of both methods is worse 
in the dynamical case which means that there might be 
more artifacts from a single method. 

Also here we find some kind of ridge line of the best 
matching values of S. Although this ridge is quite pro- 
nounced for weak filtering, we are confronted with the 
problem to identify the best matching parameters in the 
strong filtering regime, i.e. for less than 50 Laplace modes 
and more than 20 smearing steps. Furthermore, there is 
no set of best matching parameters with S > 0.5 in this 
regime and thus one can hardly say that both methods 
agree at all. 

As another measure for the similarity of topological 
landscapes at best matching we use the relative point over- 
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Figure 2: Comparison of APE with Stout smearing and APE smearing with Laplace filtering for quenched and dynamical configurations. 
The degree of similarity S (see the definition |7|) is given by contour lines. Shown are mean values of 80 configurations for APE vs. Stout 
smearing and 30 configurations for APE smearing vs. Laplace filtering. The comparisons APE vs. Stout smearing for quenched (upper left) 
and dynamical configurations (upper right) give very similar results, only the band with S > 0.95 is slightly broader in the quenched case. 
The values of S for APE smearing vs. Laplace filtering does not reach such high values and the agreement is worse in the dynamical case 
(lower row). For mild filtering we have plotted a ridge line of best matching parameters that will be used in the cluster analysis later. 



lap (RPO) [3J: 

»*B= E wE 1 - ( 9 ) 

n:£X A nX B x£Xa<JX B 

qA(x)qB{x)>Q 

Here Xa is the set of points x with topological density 
1 0*0 1 above a cut-off which is adjusted such that the 
volume fraction is the same for both methods A and B, 
voI(Xa)/VoI = voI(Xb)/VoI — f. Again, a value sab 
close to I signals a good agreement between qa and qs- 

The values of the RPOs are shown in Fig. [3j They are 
rather constant over a large range of volume fractions /. 

The RPO for APE vs. Stout smearing is quite large 
supporting the similarity of the two smearing methods. 
The RPO for APE smearing and Laplace filtering is smaller 



and comparable to its values found for quenched SU(2) 
configurations [3 . Again we find the tendency of less 
agreement in the case of dynamical configurations. 

5. Cluster Analysis of the Topological Charge Den- 
sity 

In order to obtain information about the topological 
structure of the QCD vacuum that can be compared to 
continuum models, we analyze the cluster structure of the 
topological charge density. Two lattice points belong to 
the same cluster, if they are nearest neighbors and have 
the same sign of the topological charge density. 

For such clusters in SU{2) we had found an interesting 
power law [3J. To that end one cuts the absolute value 
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Figure 3: Releative point overlap sab (see Eq.[9]l of APE and Stout 
smearing (7 iterations each) as well as for APE smearing and Laplace 
filtering (7 iterations vs. 500 Laplace modes) averaged over 10 con- 
figurations each. 
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Figure 4: Number of clusters as a function of the fraction of all points 
above a variable cut-off averaged over 30 dynamical configurations. 
The exponent of the underlying power-law is highly characteristic for 
the topological structure and directly related to the size distribution 
of topological objects. 



of the topological charge at a (variable) cut-off c7 cut and 
considers the number of clusters above this cut-off as a 
function of the fraction / of all points obeying \q{x)\ > q cu t 
w.r.t. the total number of lattice points. An example of 
this power-law for SU(3) can be found in Fig. |4j 

The exponent £ = dlog Aciust.(<7cut) /d log Appoints (<7cut) 
of this power law is highly characteristic for the topological 
structure of the QCD vacuum. Different models lead to 
different predictions. This allows for a sensitive test. 

Pure noise, for instance, yields an exponent £ = 1, 
because every point forms its own cluster. On the other 
hand, the exponent is close to zero for very smooth den- 
sities with large structures, because in such a situation a 
small change in the cut-off will add new points but not 
reveal new clusters. 

The exponent can be related to the size distribution 



d(p) ~ (P of equally charged topological objects with arbi- 
trary shape function and arbitrary dimensionality [3J. For 
instance a dilute gas of instantons would have £ = 0.64 
in the quenched case and £ = 0.66 in the dynamical case 
(N f = 2). 

A cluster analysis has one great advantage. It allows 
to reduce ambiguities coming from a single filter by taking 
only those clusters into account, which are common to 
different filters. If there is an artifact coming from one 
method, it is unlikely that this artifact will also be seen 
by the other one, such that the common structures are 
almost free of ambiguities. 

Of course, we perform this analysis after having matched 
the filtering parameters as described in the previous sec- 
tion. The exponent for clusters common to APE and 



Stout smearing can be found in Fig. 5(a) Obviously, the 
exponents of the dynamical configurations lie above the 
quenched values. The difference of the cluster exponents 
quenched vs. dynamical vanishes for stronger smearing (~ 
30 steps) and the exponents settle down to almost the 
same plateau. So we have reasons to believe that too much 
smearing destroys the impact of dynamical quarks. 

In Fig. |5(b)| the same analysis is done for Laplace fil- 
tering (only). We find a slightly bigger difference of the 
exponents of dynamical and quenched configurations. Fur- 
thermore, this difference remains, unlike for APE smear- 
ing, at every stage of filtering. At first sight the value for 
50 modes seems to be odd, but this is just an artifact of the 
incomplete reconstruction of the topological background. 
For a small number of modes, we get a very spiky struc- 
ture and we will find many new clusters, if we lower the 
cut-off. Thus, we get a higher exponent. 

To get rid of the ambiguities of APE smearing and 



Laplace filtering, we perform in Fig. 5(c) a matched clus- 
ter analysis, i.e. we consider only those clusters, which are 
common to both methods. We use the optimal set of fil- 
ter parameters according to the maximal values of S (and 
corresponding to points on the ridge line of Fig. |2| on the 
abscissa. 

Fig. |5(c) shows that the cluster exponent is larger in 
the dynamical case than in the quenched case. According 
to the considerations from above, this implies a larger ex- 
ponent (3 of the size distribution. Hence, very small topo- 
logical objects become suppressed when quarks are taken 
into account. 

In the plots we have included the values of the exponent 
in instanton gases. These have the same ordering with a 
slightly smaller difference. Generally our measured values 
of the exponent are not far off the instanton gas values 
(which was not the case in SU(2) [3J). 

Another interesting quantity is the absolute number of 
clusters above a certain cut-off, see Fig. [6j In the dynam- 
ical case we find almost twice as many clusters. Hence, 
we observe a more fragmented topological structure in the 
presence of dynamical quarks. 

Our findings can be compared to the results obtained 
by the Adelaide group 0. They have observed an "in- 
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Figure 5: Exponent £ for the clusters found by APE smearing (a), 
Laplace filtering (b) and for a matched filtering of both methods 
(c). The exponents predicted by the dilute instanton gas have been 
included in all plots. In the latter plot the measured exponents have 
a step-like behavior, because the number of smearing steps is an 
integer quantity and we have to match several numbers of Laplace 
modes to the same number of smearing steps. The smooth lines are 
interpolations corresponding to a "continuous" smearing parameter. 
The number of smearing steps (in parentheses) refers to the matching 
of quenched configurations and is slightly different for the dynamical 
case. Errors from the ensemble average have been included, but are 
partly too small to be seen. 
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Figure 6: Total number of distinct clusters for a constant fraction / = 
0.0755 of points lying above the cut-off with smooth interpolation. 

creasing density of nontrivial field configurations" and a 
suppression of small instantons. Using a completely dif- 
ferent approach, we come to the same conclusion. The ad- 
vantages of our analysis are that our dynamical fermions 
have nicer chiral properties than the staggered ones. More- 
over, we do not have to postulate any shape function of the 
topological objects. Therefore, our method works for all 
topological building blocks and with our matched analysis 
we can obtain results which are almost free from ambigu- 
ities of the filtering process. 

6. Conclusions 

We have performed a comparative study on the local 
topological structure of quenched and dynamical config- 
urations seen by various filtering methods, namely APE 
smearing, Stout smearing and Laplace filtering. 

We found an almost one-to-one correspondence of the 
topological charge densities resulting from APE and Stout 
smearing, if the filtering parameters are chosen properly. 
This extends findings in lattice perturbation theory [8] and 
on non-topological observables for few smearing steps [5]. 
Our results reflect the local similarity of both methods for 
up to 50 steps and their strongly correlated topological 
charge densities. 

As a consequence APE and Stout smearing can be seen 
as equal for the investigation of topological objects on the 
lattice. On the other hand, both methods will therefore 
suffer from almost the same artifacts and other, indepen- 
dent methods are needed to get results free from ambigu- 
ities. 

To this end Laplace filtering, which is based on a trun- 
cated spectral expansion of the Laplacian, has been in- 
vestigated in addition to smearing. In the case of mild 
filtering, we found a qualitative agreement of Laplace fil- 
tering and APE smearing. This is in accordance with the 
quenched SU{2) results of Bruckmann et al. [3]. For 
SU(3) the agreement is slightly better in the quenched 
case. However, for the strong filtering case we have found 
a clear disagreement, which is even more pronounced in 
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the dynamical case. This finding flags a warning that any 
kind of smearing should be applied with great care. 

In the central part of this work the topological clusters, 
revealed by the different methods, have been investigated. 
To this end we have analyzed the power-law behavior of 
topological charge clusters seen by the individual methods 
and common to both methods. 

This analysis shows clearly a larger exponent in the 
presence of dynamical quarks. In order to interpret this 
observation, we can relate the cluster exponent to the ex- 
ponent of the size distribution d{p) ~ p* 3 of topological ob- 
jects in corresponding models. Our findings give a larger 
coefficient j3 in the dynamical case, hence, very small topo- 
logical objects become suppressed. 

Moreover, the total number of clusters for a constant 
total cluster volume is substantially higher for dynamical 
configurations. Thus the topological structure of the vac- 
uum is much more fragmented in the presence of fermions. 
It would be interesting to further investigate the nature of 
these topological lumps. 
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